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Abstract 

In this paper the acoustic field propagating in the early hot {p = e/3) universe of arbitrary space 
curvature {K = 0, ±1) is considered. The field equations are reduced to the d'Alembert equation 
in an auxiliary static Roberson- Walker space-time. Symbolic computation in Mathematica is 
applied. 
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1. Introduction 

The Sachs- Wolfe theorem (Sachs and Wolfe, 1967, pp. 76-77) contains two separate 
theses formulated for two different equations of state: the first ((i) p. 76) for the pres- 
sureless matter {p = 0), and the second ((ii) p. 77) for ultrarelativistic gas (p — e/3). 
The first thesis was recently recalculated and rediscussed by Sachs et al., (2007). In this 
paper we will concentrate on the second, and for the purposees of this paper we will call 
it the acoustic theorem to clearly distinguish between them both. 

Acoustic theorem refers to the spatially flat (K = 0), hot {p = e/3) Friedman- 
Robertson- Walker universe, and to the scalar perturbation propagating in it. It satates, 
that by the appropriate choice of the perturbation variable one can express the prop- 
agation equation in the form of d'Alembert equation in Minkowski space-time. Scalar 
perturbations in the flat, early universe propagate in a similar manner to electromag- 
netic or gravitational waves (Sachs and Wolfe, 1967, p. 79, see also Grishchuk 2007). 
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On the other hand, the d'Alembert equation describing scalar field in the dust {p = 
0) cosmological model, can be transformed into the d'Alembert equation in the static 
Robertson- Walker space- time regardless of the universe space curvature (Klainerman and 
Sarnak, 1981). We may, therefore, suppose that the flatness assumption in the Sachs- 
Wolfe theorem is excessive and that the theorem is true in the general case. 

Indeed, as we show below, the scalar perturbations reveal similar features in both, the 
flat and curved universes, i.e. they propagate as waves in the static Robertson-Walker 
space-time. The proof of this fact is the subject of this paper. The Sachs Wolfe acoustic 
theorem may apply to the microwave background radiation (Sachs et al., 2007). 

The paper consists of a brief description of the problem, and the short code in Math- 
ematica. Throughout this paper c = 1 and SttG = 1. 



2. Scalar perturbations on the Friedman Robertson Walker background 

We adopt Robertson-Walker metrics in spherical coordinates x'^ = [r], x, <p} 
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with the scale factor a(ry) appropriate for the equation of state p = e/3 

sin ^y/Kri^ 



a{ri) = 
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Consider the scalar perturbations in the synchronous system of reference. The metric 
tensor correction is then determined by two scalar functions £ and C (Stewart, 1990) 



Sg^J.Q = 0, 



(3) 
(4) 



A stands for the Beltrami-Laplace operator on the r] = const hypersurface. Eventually, 
the metrics takes the form 

Qtiv = 9{nw)iiv + Sgiiv (5) 
In the first order of the perturbation expansion the Einstein equations with the metric 
tensor (5) reduce to 



drj 
52 
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where we put /i = C, A = — A£ (to fit to (Lifshitz and Khalatnikov, 1963, Landau and 
Lifshitz, 2000)). 0(77) stands for sound velocity. The energy density contrast S = Se/e 
reads 

1 



3e(ry)a2(7?) 



i^^fiix^ - {3K + A) [A(x") + fi{x'')] 



(8) 



a{ri) drj 

Let us define a new perturbation variable ^ with help of the second order differential 
transformation of the density contrast 5 
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The function 'F(x°') is the solution of the d'Alcnibert equation 

|^^(x'^)-1ai^(.t-) = 
with the Beltrami-Laplace operator A — ^'^•'(7mnV™V" acting in the space 
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Theorem 1 (Sachs- Wolfe acoustic theorem for K = 0,±1). Scalar perturbations in 
the hot (p — tji) Friedman-Robertson-Walker universe of arbitrary space curvature 
(K = 0, ±1 ), expressed in terms of the perturbation variable ^ (9) obey the wave equation 
(10) in the static Robertson-Walker space-time g = diag{~l, ^^■'(7). 



Proof. 

We perform this calculation in Mathematica. The code is enclosed below, and also can 
be downloaded from: 



(* http://drac.oa.uj.edu.pl/usr/woszcz/kody/acoucticRWcode.nb *) 



(* Coordinates *) 
Clear[X , o, c, e, H, Lap, 5*] 
X = {77, X, 7?, XX = {rj_, X-, ^-}; 
x = Sequence@@X; xx = Sequence@@XX; 

(* Friedman-Robertson- Walker background for e = 3p *) 



c[rj_ 



] = \ Vi-\ = Limit 



K 



H[V-] - ^//Simplify; e[,7.] = -3^ + SiJ^ V/Simplify; 

(* Lifshitz-Khalatnikov perturbation equations *) 
^(2,0,0,0) [xx] = -^a„A[x] - i(Lap[A][x] +Lap[Ai][x]); 
^(2,0,0,0) [xx] = - (2 + ScM^) ^ d,fr[x] + i (l + Sc^^) 



( (Lap [A] [x]+3KX[x]) + (Lap ^ [x] + 3X^ [a;] ) ) ; 
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(* The density contrast *) 
-^^=3™^ (-((Lap[A] W + 3KX[x]) + (Lap[M]H + Si^^N)) + 3^5{,,i}/xN) 



//Simplify; 



(* Higher order derivatives *) 

;^(kJ„teger/;fe>2,0,0,0)[^^].^^ [;^(2,0,0,0) [^] ^ {^^ ^ _ 2}] 
^(k.I„teger/;fc>2,0,0,0) [^x] :=D [/X^^.O.O'O) [x] ,{rj,k-2}] 

^(2,di-,d2-,d3.)[xx] = D [A(2.o.o.o)[x],{x,dl},{i?,d2},{<^,d3}] ; 
^(2,di-,d2.,d3.)[^j = D [/x(2.0A0)[^]^|^^dl},{^?,d2},{^,d3}] ; 

(* Laplacian in the maximalfy symmetric 3-dim curved space *) 
w[K_, X-] = Limit [2yKCot [-/kx] ^I"^ ^ K]; 

Lap[L][xx] = j^p^ (Csc[i?]2a{^,2}/N +Cot[,9]9„/[x] +a{^,2}/N) 

+ w[K,x]d^f[x] +9{^,2}/N; 

(* The gauge-invariant variable ^ *) 



r[K_, = Limit 



(* The wave equation *) 
5{^,2}*W - iLap[4'][a;]==0; 

Timing[Simplify[%,TimeConstraint 5000]] 

□ 



3. Discussion 

By obtaining equation (10) we have expressed the problem of cosmological density 
perturbations in terms of the semiclassical theory (the field theory in curved space- 
time — sec BirrcU and Davies, 1982). Consequently, the key facts resulting from the 
semiclassical theory apply also to the perturbed expanding universe. 

Scalar perturbations in the early universe form the acoustic field. Acoustic waves are 
dispersed by the space curvature. The dispersion relation for the equation (10) takes the 
form (Birrell and Davies, 1982, eq. (5.27)) 



and the group velocity reads 



lik'-K) (12) 



d k 

dk'^^Tw^^- ^^^^ 



In flat space (K = 0) the group velocity is constant and equal to l/VS. This fact has 
been already established by Sachs and Wolfe (Sachs and Wolfe, 1967). 
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Beyond the K = case the dispersion relation is nonhnear. In the space of negative 
curvature K = —1 the waves behave as the scalar field with the mass m = 1. The group 
velocity is the function of k and decreases to zero in the k ^ limit 

hm Wg = (14) 

while the limit frequency is still positive 

lim uj ~ —=. (15) 
fc^o V3 

Acoustic behaviour does not extend to the solutions of imaginary k: k^ G { — 1,0} (sup- 
plementary series (Gelfand et al., 1966, Neimark, 1964), supercurvature modes (Lyth and 
Woszczyna, 1995)). 

In the closed universe (K = 1) the wave numbers are integer and 1 < fc . The general 
solution is a countable combination of the hyperspherical functions. 

On basis of ^{x'^) one can reconstruct the solution for the density contrast 5{x'') 
by means of the reciprocal relation (9). (Compare analogue procedure for scalar field 
(Klainerman and Sarnak, 1981). By reconstructing 5{x'^) one is forced to restore the 
undesirable phenomenon of gauge ambiguity. 
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